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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement. 

(Symbols  of  SI  units  given  in  parentheses) 


To  convert  from 

to 

Multiply  by 

angstrom 

meters  (m) 

1.000  000  X 

atmosphere  (normal) 

kilo  pascal  (kPa) 

1.013  25  X  E 

bar 

kilo  pascal  (kPa) 

1.000  000  X 

barn 

p  2 

meter  (m  ) 

1.000  000  X 

British  thermal  unit 
(thermochemical ) 

joule  (J) 

1.054  350  X 

calorie  (thermochemical) 

joule  (J) 

4.184  000 

P 

cal  ( thermochemi cal )/cm 

2  2 
mega  joule/m  (MJ/m  ) 

4.184  000  X 

curie 

giga  becquerel  (GBq)* 

3.700  000  X 

degree  (angle) 

radian  (rad) 

1.745  329  X 

degree  Fahrenheit 

degree  kelvin  (K) 

=  (t°  f  +  459. 

electron  volt 

joule  (J) 

1.602  19  X  [ 

erg 

joule  (J) 

1.000  000  X 

erg/second 

watt  (W) 

1.000  000  X 

foot 

meter  (m) 

3.048  000  X 

foot-pound-force 

joule  (J) 

1.355  818 

gallon  (U.S.  liquid) 

meter^  (m^) 

3.785  412  X 

inch 

meter  (m) 

2.540  000  X 

jerk 

joule  (J) 

1.000  000  X 

joule/kilogram  (J/kg)  (radiation 
dose  absorbed) 

Gray  (Gy)** 

1.000  000 

kilotons 

terajoules 

4.183 

kip  (1000  Ibf) 

newton  (N) 

4.448  222  X 

P 

kip/inch  (ksi) 

kilo  pascal  (kPa) 

6.894  757  X 

ktap 

2  2 
newton-second/m  (N-s/m  ) 

1.000  000  X 

micron 

meter  (m) 

1.000  000  X 

mi  1 

meter  (m) 

2.540  000  X 

mile  (international) 

meter  (m) 

1.609  344  X 

ounce 

kilogram  (kg) 

2.834  952  X 

pound-force  (Ibf  avoirdupois) 

newton  (N) 

4.448  222 

pound-force  inch 

newton-meter  (N*m) 

1.129  848  X 

pound- force/inch 

newton/meter  (N/m) 

1.751  268  X 

pound-force/foot 

kilo  pascal  (kPa) 

4.788  026  X 

E  -10 
+2 

E  +2 
E  -28 

E  +3 

E  -2 
E  +1 
E  -2 
67)71.8 
-19 
E  -7 
E  -7 
E  -1 

E  -3 
E  -2 
E  +9 


E  +3 
E  +3 
E  +2 
E  -6 
E  -5 
E  +3 
E  -2 

E  -1 
E  +2 
E  -2 
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Conversion  factors  for  US.  customary  to 
metric  (SI)  units  of  measurement.  (Continued) 


(Symbols  of  SI  units  given  in  parentheses) 


To  convert  from 

to 

Multinlv  hv 

pound -force/inch^  (psi) 

kilo  pascal  (kPa) 

— _ 

6.894 

- ZaU 

757 

pound-mass  (ibm  avoirdupois) 

p 

kilogram  (kg) 

4.535 

924  X 

E  -1 

pound-mass-foot 

(moment  of  inertia) 

p  p 

kilogram-meter  (kg -nr) 

4.214 

on  X 

E  -2 

pound-mass/foot 

kilogram/meter^  (kg/m^) 

1.601 

846  X 

E  +1 

rad  (radiation  dose  absorbed) 

Gray  (Gy)** 

1.000 

000  X 

E  -2 

roentgen 

coulomb/kilogram  (C/kg) 

2.579 

760  X 

E  -4 

shake 

second  (s) 

1.000 

000  X 

E  -8 

slug 

kilogram  (kg) 

1.459 

390  X 

E  +1 

torr  (mm  Hg,  0°  C) 

kilo  pascal  (kPa) 

1.333 

22  X 

E  -1 

*The  becquerel  (Bq)  is  the  SI 

unit  of  radioactivity;  1  Bq  = 

1  event/s. 

**The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 

A  more  complete  listing  of  conversions  may  be  found  in  "Metric  Practice  Guide 
E  380-74,"  American  Society  for  Testing  and  Materials. 
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SECTION  1 
INTRODUCTION 

The  purpose  of  the  Cylindrical  In-Situ  Test  (CIST),  developed  in  1971,  is 
to  measure  the  dynamic  response  of  geologic  materials  to  a  cylindrically 
symmetric  high  explosive  shock  input.  This  in-situ  material  properties  test 
was  developed  at  the  Air  Force  Weapons  Laboratory  (AFWL)  to  eliminate  some  of 
the  shortcomings  inherent  in  the  material  sampling/laboratory  testing  approach 
to  material  property  determination.  A  vertical  cylindrical  load  geometry  is 
used  because  of  its  unique  capability  for  simultaneous  loading  of  a  number  of 
near-surface  horizontal  layers.  Resulting  data  are  then  used  to  determine  the 
in-situ  material  properties  of  each  layer. 

In  the  past  few  years,  extensive  theoretical  and  experimental  programs 
have  been  undertaken  in  an  effort  to  quantify  the  response  of  strategic 
structures,  sited  in  widely  varying  geologies,  to  various  airblast  and  ground 
shock  loadings.  Characterization  of  this  response  depends  on  an  adequate 
mathematical  model  of  the  material  surrounding  the  structure;  therefore,  a  prime 
requirement  for  determining  the  response  of  strategic  structures  to  nuclear  and 
high  explosives  blast  loadings  is  an  understanding  of  the  response  of  the 
surrounding  geologic  materials  to  impulsive-type  loads. 

In  order  to  obtain  this  dynamic  response  information,  the  CIST  18 
experiment  was  detonated  on  the  HAVE  FIOST  Test  Site,  located  on  Luke  Air  Force 
Range,  on  22  October  1976  (HAVE  HOST  being  the  name  of  the  AFWL/SAMSO/DNA  test 
series  to  validate  MX  protective  structure  concepts).  The  material  models 
derived  from  the  CIST  data,  along  with  laboratory  data,  will  be  used  in  the 
structural  response  and  ground  shock  calculations  for  the  HAVE  HOST  Test  Series, 
and  in  the  more  general  Multiple  Aim  Point  ground  shock  prediction  calculations. 

In  addition,  the  CIST  18  data  yi el ds  important  i nformation  concerning 
Multiple  Aim  Point  geologies.  Previous  CISTs  have  been  conducted  on  the  Nevada 
Test  Site  (CISTs  5-7)  (Ref.  1)  and  White  Sands  Missile  Range  (CISTs  15  and  16) 
(Ref.  2).  These  six  CIST  events  provide  valuable  information  on  a  wide  range 
of  materials  that  occur  in  the  southwestern  United  States,  ranging  from 
saturated  clays  to  dry  sands  to  hard  rock. 
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SECTION  2 
TEST  DESCRIPTION 

The  General  Test  Plan  for  the  Cylindrical  In-Situ  Test  (CIST),  AFWL-TR- 
74-136  (Ref.  3),  provides  the  general  information  on  the  fielding  of  CIST  tests 
and,  unless  otherwise  noted,  no  procedural  deviations  were  made  for  CIST  18. 

CIST  18  was  located  on  the  eastern  boundary  of  the  HAVE  HOST  Test  Site 
(Figure  2.1),  approximately  30  miles  southeast  of  Yuma,  Arizona,  on  the  Luke 
Air  Force  Range.  Figure  2.2  is  a  map  of  Arizona,  giving  the  relative  locations 
of  Yuma  and  the  HAVE  HOST  Test  Site. 

The  major  difference  between  the  HAVE  HOST  CIST  and  p’^_;ious  CISTs  is  that 
two  separate  events,  18D  and  18S,  were  detonated.  The  two  CISTs  were 
necessitated  by  the  inability  to  find  a  single  location  that  contained  all  the 
materials  needed  to  characterize  the  HAVE  HOST  Test  Site.  These  materials  range 
from  poorly  to  moderately  cemented  alluvium  and  highly  jointed  and  fractured 
granite  to  highly  jointed  competent  granite. 

CIST  18D  (deep),  for  which  geology  and  instrumentation  plan  is  given  in 
Figure  2.3,  was  the  original  site,  but  inspection  of  the  undisturbed  samples 
indicated  a  lack  of  the  poorly  cemented  surface  material  prevalent  over  much  of 
the  HAVE  HOST  Test  Site.  CIST  18S  (shallow),  shown  in  Figure  2.4,  was  located 
150  feet  northwest  of  18D  in  a  drainage  channel  where  an  abundance  of  poorly 
cemented  alluvial  material  existed.  A  profile  of  the  subsurface,  interpreted 
from  a  refraction  survey,  on  a  line  between  the  two  CIST  ground  zeroes  (GZ)  is 
shown  in  Figure  2.5. 

2.1  SITE  GEOLOGY 

The  CIST  18  experiment  was  originally  designed  to  test  two  different 
materials--a  hard,  jointed  granite  bedrock  and  a  dry  alluvium  formed  from 
the  mechanical  decomposition  of  the  granite.  Inspection  of  the  undisturbed 
samples  at  the  18D  GZ  revealed  three  different  types  of  silty  sandy  alluvium  and 
two  different  types  of  granite. 

The  silty  sand  is  present  over  the  granite  throughout  the  entire  HAVE  HOST 
Test  Site,  except  where  the  granite  outcrops.  The  only  apparent  physical 
difference  in  the  alluvium  is  the  amount  of  cementation  present,  and  there  is  no 
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Figure  2.1.  HAVE  HOST  test  site. 
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Figure  2.3.  CIST  18D  generalized  geology  and  instrumentation  plan. 
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Figure  2.5.  Geology  cross  section  between  CIST  GZs. 


AFWL-TR-77-81 


trend  of  increasing  cementation  with  depth.  There  is  some  evidence  that  the 
degree  of  cementation  is  reduced  as  the  grain  size  of  the  alluvium  inci'eases. 

The  alluvium  was  logged  as  "sand,  fine  to  coarse,  silty  with  small  pebbles." 
The  pebbles  were  almost  exclusively  weathered  granite.  Once  the  alluvium  had 
been  disturbed,  it  was  impossible  to  visually  determine  the  depth  or  the  location 
at  which  the  sample  originated.  In  general,  three  materials  were  found  to  exist; 
a  slightly  to  moderately,  moderately,  and  moderately  to  well  cemented  alluvium. 

For  the  purposes  of  this  report,  slightly  cemented  is  defined  as  easily  broken 
apart  with  the  hands.  Moderately  cemented  is  defined  as  broken  apart  with  the 
hands  with  some  difficulty,  and  weTl  cemented  is  defined  as  impossible  to  break 
apart  with  hands  or  broken  apart  with  hands  with  great  difficulty.  One  of  the 
major  purposes  of  CIST  18  was  to  determine  the  effect  of  these  varying  degrees 
of  cementation  on  the  material  response. 

In  the  180  explosive  hole,  the  soil  extended  to  the  granite  bedrock  at  a 
depth  of  53.9  feet,  and  the  transiti^on  zone  between  the  alluvium  and  the  granite 
was  very  pronounced  and  easily  identifiable.  The  granite  was  competent  and 
highly  jointed  with  a  prominent  north-south  joint  trend.  The  jointing  was 
evident  in  the  granite  outcrop  100  meters  east  of  the  test  site  where  the  joint 
spacing  was  less  then  3  inches.  This  close  joint  spacing  caused  difficulties  in 
obtaining  undisturbed  rock  cores.  At  certain  depths  in  the  granite,  the  method 
of  coring  (Longyear  NQ  wireline  core  barrel)  pulverized  the  material  so. badly 
that  samples  longer  that  1-1/2  inches  were  rare. 

Inspection  of  the  intact  cores  revealed  numerous  joints  that  were  healed 
well  enough  so  as  not  to  separate  on  sampling.  The  difference  in  obtaining  intact 
and  fractured  cores  appeared  to  be  the  strength  of  the  joint  cementation.  This 
observation  also  leads  to  the  two  different  types  of  granite  tested--highly 
fractured  weakly  cemented,  and  less-fractured  stronger-cemented.  As  with  the 
alluvium,  no  trend  of  increasing  competency  with  depth  was  found,  although  this 
result  is  to  be  expected. 

2.2  CIST  18D 

The  18D  er.plosive  cavity  was  drilled  to  a  depth  of  73.7  feet.  The  hole  was 
cased  with  a  24  inch  steel  culvert  down  to  the  alluvium/rock  interface  at  53.9 
feet.  A  22  inch  explosive  cavity  was  then  drilled  to  73.7  feet.  The  charge 
density  of  5  pounds  of  detonating  cord/linear  foot  of  cavity  remained  constant 
in  both  the  rock  and  alluvium. 
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Except  for  a  few  feet  of  surface  material,  the  entire  soil  profile  down  to 
rock  was  moderately  well  cemented.  The  gauges  at  the  15.0  feet  level  were  placed 
in  moderately  to  well  cemented  alluvium,  and  the  gauges  at  the  31.0  feet  level 
were  placed  in  moderately  cemented  alluvium.  The  gauges  at  the  62.0  feet  level 
were  in  the  most  competent  material  available  (highly  jointed,  competent  granite). 

2.3  CIST  18S 

The  shallow  explosive  cavity  was  20.0  feet  deep  and  lined  with  a  24  inch 
steel  culvert.  Except  for  a  1/2-foot  hard  crust,  the  first  3  feet  of  material  were 
cohesionless.  When  drilling  instrumentation  holes  through  this  material,  bulbous 
cavities  3  to  4  feet  in  diameter  were  formed  under  the  surface  crust.  This 
necessitated  casing  the  first  4  to  5  feet  to  prevent  fall -in.  After  the 
explosive  cavity  was  augered,  the  24  inch  steel  culvert  was  grouted  in.  The 
materials  below  this  loose  surface  deposit  were  noticeably  less  cemented  than  in 
18D  and  were  logged  as  slightly  to  moderately  cemented. 

2.4  INSTRUMENTATION 

The  standard  instrumentation  packages,  accelerometer  canisters  and  stress 
gauge  paddles  were  used  in  both  18D  and  18S.  However,  placement  techniques  were 
modified  to  be  compatible  with  the  geologic  materials  to  be  tested. 

With  the  exception  of  the  18D  undisturbed  rock  cores  and  the  section  of  the 
explosive  cavity  that  was  in  granite,  all  the  holes  were  drilled  dry.  The 
instrument  holes  were  bored  using  compressed  air  as  the  drilling  fluid,  and  the 
part  of  the  explosive  cavities  in  alluvium  were  excavated  using  a  bucket  auger. 
This  was  done  to  prevent  water  intrusion  into  the  walls  of  the  hole  which  might 
alter  normal  in-si tu  moisture  contents  and,  therefore,  seriously  affect  the  data. 

The  procedure  used  in  the  excavation  of  the  explosive  cavity  was  as  follows: 
A  27  inch  bucket  auger  was  used  to  remove  the  alluvial  material  down  to  the 
granite  bedrock.  A  24  inch  steel  casing  was  then  placed  in  the  cavity  and 
grouted  in.  After  the  grout  had  hardened,  the  section  of  the  explosive  cavity 
in  the  rock  was  drilled  using  a  22  inch  rock  bit  and  drilling  mud.  The  22  inch 
bit  was  chosen  as  the  largest  bit  that  could  be  fit  inside  the  culvert. 

Figures  2.6  and  2.7  give  the  plan  views  of  the  18S  and  18D  instrument  holes. 
Table  1  gives  the  "As-Built"  measurement  list. 
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Table  2.1  CIST  18  Measurement  List 


Measurement 

Number 

Measurement  Designation 

Gauge  Range 

Cal ibration 

01 

CIST-18S-F-E-L13- 12. 0-45-3. 0-A-V 

10,000 

g's 

2273.0 

g's 

02 

CIST-18S-F-E-L13-12. 0-45-3. 0-A-H 

10,000 

g's 

4936.9 

g's 

03 

CIST-18S-F-E-L20-2. 5-25-5. 0-SE-H 

2,000 

psi 

2934. C 

psi 

04 

CIST-18S-F-E-L14-12. 0-90-5. 0-A-V 

2,000 

g‘s 

697.2 

g's 

05 

CIST-18S-F-E-L14-12. 0-90-5. 0-A-H 

5,000 

g's 

1582.4 

g's 

06 

CIST-18S-F-E-L17-12. 0-270-5. 0-SE-H 

2,000 

psi 

1344.6 

psi 

07 

CIST-18S-F-E-L15-12. 0-135-8. 0-A-V 

2,000 

183.4 

g's 

08 

CIST- 18S-F-E-L15-12. 0-135-8. 0-A-H 

2,000 

g's 

176.5 

g's 

09 

CIST-18S-F-E-L19-12. 0-0-8. 0-SE-H 

2,000 

psi 

384.6 

psi 

10 

CIST-18D-F-E-L02-15. 0-330-3. 0-A-V 

10,000 

g's 

2778.0 

g's 

11 

CIST-18D-F-E-L02- 15. 0-330-3. 0-A-H 

10,000 

g's 

5014.8 

g's 

12 

CIST- 18D-F-E-L09-15. 0-60-3. 0-SE-H 

2,000 

psi 

3046.6 

psi 

13 

CIST-18D-F-E-L03-15. 0-0-5. 0-A-V 

2,000 

g's 

711.8 

g's 

14 

CIST-18D-F-E-L03-15. 0-0-5. 0-A-H 

5,000 

g's 

1642.5 

g's 

15 

CIST-18D-F-E-L10-15. 0-90-5. 0-SE-H 

2,000 

psi 

1290.6 

psi 

16 

CIST-18D-F-E-L04-15. 0-30-8. 0-A-V 

2,000 

g's 

180.2 

g's 

17 

CIST-18D-F-E-L04-15, 0-30-8. 0-A-H 

2,000 

g's 

178.4 

g's 

18 

CIST-18D-F-E-L11-15. 0-120-8. 0-SE-H 

2,000 

psi 

408.9 

psi 

19 

CIST- 18S-F-E-L16-12. 0-225-5. 0-SE-T 

2,000 

psi 

1259.1 

psi 

20 

CIST-18S-F-E-L18-12. 0-315-8. 0-SE-H 

2,000 

psi 

3006.3 

psi 

21 

CIST-18D-F-E-L02-31. 0-330-3. 0-A-H 

10,000 

g's 

4954.2 

g's 

22 

CIST-18D-F-E-L05-62. 0-150-3. 0-A-H 

30,000 

g's 

25259.6 

g's 

23 

.  CIST-18D-F-E-L05-69. 0-150-3. 0-A-V 

30,000 

g's 

11630.0 

g's 

24 

CIST-18D-F-E-L05-69. 0-150-3. 0-A-H 

30,000 

g's 

31111.1 

g's 

25 

CIST-18D-F-E-L03-31. 0-0-5. 0-A-H 

5,000 

g's 

1681.0 

g's 

26 

CIST-18D-F-E-L06-62. 0-210-5. 0-A-H 

10,000 

g's 

4822.0 

g's 

27 

CIST- 18D-F-E-L06-69. 0-210-5. 0-A-V 

5,000 

g's 

2555.2 

g's 

28 

CIST-18D-F-E-L06-69. 0-210-5. 0-A-H 

20,000 

g's 

10330.0 

g's 

29 

CIST- 18D-F-E-L04-31. 0-30-8. 0-A-H 

2,000 

g's 

177.4 

g's 

30 

CIST- 18D- F-E-L07-62 . 0-240-8. 0-A-H 

5,000 

g's 

1260.0 

g's 

31 

CIST-18D-F-E-L07-69. 0-240-8. 0-A-V 

2,000 

g's 

674.8 

g's 
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Table  2.1  CIST  18  Measurement  List  (Continued) 


Measurement 


Number 

Measurement  Designation 

Gauge  Range 

Cal i bration 

32 

CIST-18D-F-E'L07-69.0-240-8.0-A-H 

5,000  g's 

2458.3 

g's 

33 

CIST-18D-F-E-L08-62. 0-270-12. 0-A-H 

2,000  g's 

341.3 

g's 

34 

CIST-18D-F-E-L08-69. 0-270- 12. 0-A-V 

2,000  g's 

209.7 

g's 

35 

CIST- 18D-F-E-L08-69. 0-270-12. 0-A-H 

2,000  g's 

698.4 

g's 

36 

CIST-18D-F-E-L12-12.0-n-1.0-CP-V 

10,000  psi 

9813.9 

psi 

All  channels  redundantly  recorded,  except  for  MN  10,16. 
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SECTION  3 
DATA  AND  ANALYSIS 

Figure  3.1  is  a  plot  of  peak  horizontal  particle  acceleration  versus  range 
for  all  the  CIST  18  alluvium  layers  (both  18S  and  18D).  From  this  plot,  the 
expected  material  response  due  to  cementation  effects  is  not  obvious.  These 
values  plot  within  the  scatter  of  the  CIST  5  (Area  10,  NTS)  alluvial  data.  The 
predictions  for  the  alluvial  response  were  based  on  the  CIST  16  (Queen  15  Area 
WSMR)  and  were  adequate. 

A  plot  of  maximum  horizontal  velocity  versus  range  (Figure  3.2)  shows 
virtually  no  difference  in  the  peak  velocities  of  the  three  layers.  This  result 
reinforces  che  previously  stated  observation:  Cementation  effects  had  a 
negligible  effect  on  maximum  values. 

Figure  3.3  is  a  plot  of  all  the  dry,  sandy  alluvium  CIST  data.  These  data 
consist  of  CIST  5  (silty  sands),  CIST  16  (clayey,  silty  sands)  and  CIST  18  (silty 
sands).  It  is  seen  that  the  CISTs  5  and  18  peak  velocity  response  is  about  the 
same,  and  the  CIST  16  dry  alluvium  data  can  be  used  as  an  upper  bound  of  all  MAP 
dry,  sandy  alluvium  data  collected  to  date. 

An  evaluation  of  the  "one-dimensionality"  of  the  alluvium  data  is  accompli  shed 
in  a  table  of  the  ratio  of  peak  vertical  particle  velocity  to  peak  horizontal 
particle  velocity  for  all  depths  and  ranges  in  the  CIST  18  experiment. 

3'  5'  8' 

18D  --15'  depth  0.093  0.123  0.018 

18S  -  12'  depth  0.023  0.049  0.026 

Because  ratios  are  on  the  order  of  0.1  or  less,  we  feel  that  the  material 

model  development  can  be  treated  as  a  one-dimensional  problem. 

The  actual  rock  response  was  about  an  order  of  magnitude  below  the  pretest 
prediction,  as  shown  in  a  plot  of  peak  horizontal  acceleration  versus  range  for 
the  jointed  granite  (Figure  3.4).  This  over-prediction  is  attributed  to  two 
factors;  the  conservative  method  of  prediction  and  an  underestimation  of  the  in- 
situ  jointing  and  fracture  conditions  of  the  granite.  The  conservative  method 
of  prediction  was  dictated  because  there  were  no  hard  rock  CIST  data  prior  to 
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Figure  3.2.  Plot  of  peak  horizontal  particle  velocity 
versus  range  for  CIST  18  alluvium. 
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Figure  3.3.  Plot  of  peak  horizontal  particle  velocity 
versus  range  for  MX  sandy  dry  alluvium. 
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CIST  18,  and  such  data  were  a  primary  goal  of  CIST  18.  Consequently,  to  insure 
that  data  were  obtained,  the  predictions  were  deliberately  conservative.  In 
addition,  two  different  sets  of  gauge  ranges  were  used  in  the  two  different 
rock  layers  and  redundant  recordings  of  all  measurements  were  made.  The  g>anite 
response  predictions  were  based  on  the  DATEX  II  granite  (Ref.  4)  (Cedar  City 
tonalite).  The  seismic  velocity  of  this  material  is  11,000  feet/second,  as 
compared  to  the  9,300  feet/second  HAVE  HOST  oranite.  The  matrix  structure  of 
the  DATEX  II  granite  is  weaker  vhat  HAVE  HOST;  however,  the  DATEX  II  joint 
spacing  is  on  the  order  of  1  to  3  feet,  whereas,  the  HAVE  HOST  joint  spacing  is 
about  3  inches.  In  addition,  the  extensive  fracturing  that  is  common  at  HAVE 
HOST  is  absent  at  the  DATEX  test  site,  making  the  tonalite,  in  a  macroscopic 
sense,  much  more  competent  than  the  HAVE  HOST  granite.  The  jointing  and 
fracturing  of  the  rock  would  tend  to  fiUer  out  high  frequency  response,  as  the 
CIST  18  granite  data  illustrate,  and  greatly  increase  the  spatial  attenuation 
that  would  occur.  In  addition,  DATEX  II  was  a  plane  wave  simulation,  and  CIST 
tests  are  a  cylindrical  wave  simulation  iiiaking  scaling  from  one  to  the  other 
extremely  difficult. 

A  plot  of  soil  stress  data  versus  range  (Figure  3.5)  indicates  the  anomalous 
nature  of  the  soil  stress  measurements.  Some  of  the  gauges  appear  to  have  over 
registered.  Possible  causes  are  faulty  placement  techniques  or  the  soil  grain- 
gauge  contact.  With  sandy  materials,  a  nonuniform  soil  contact  across  the  gauge 
face  is  possible.  If  this  happens,  very  high  point  stresses  could  develop  and 
make  the  gauge  reading  erroneous.  As  further  evidence,  the  pCpV  calculations 
assuming  a  density  (p)  of  115  pounds/feet^  and  loading  wave  speed  (Cp)  =  809 
feet/second  are  a  lower  bound  to  the  stress  data.  There  were  also  instrumenta¬ 
tion  difficulties  that  have  not  been  resolved,  i.e.,  calibration  levels  of 
several  of  the  18S  gauges  that  cast  further  doubt  on  the  validity  of  the 
measurements.  These  problems  have  prompted  the  fielding  of  a  Placement  Evalua¬ 
tion  Technique  (PET)  test  in  February  1977  at  the  HAVE  HOST  Test  Site.  Until 
the  calibration  problems  are  resolved  and  the  PET  test  data  analyzed,  the  stress 
data  should  be  considered  unusable  or  highly  suspect. 

The  time  of  arrival  (TOA)  versus  range  are  displayed  in  Figure  3.6.  These 
data  were  computed  from  the  first  motion  of  the  accelerometers.  The  seismic  ve¬ 
locities  of  the  alluvium  layers  are  2000,  3200  and  3500  feet/second  for  18D,  15 
feet;  18S,  12  feet;  and  180,  31  feet,  respectively.  It  seems  that  the  initial 
assumption  of  18S  being  less  cemented,  is  not  supported  by  the  TOA  data.  The 
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seismic  velocity  is  a  possible  indication  of  the  alluvial  cementation,  and  it  now 
appears  logical  to  assume  that  the  present  surface  water  drainage  channels  at  the 
HAVE  HOST  Test  Site  are  not  necessarily  in  the  same  position  as  they  were  in 
earlier  geologic  times  when  the  subsurface  materials  were  deposited. 

It  appears  that  soil  cementation  is  a  function  of  two  distinct  parameters; 
in-situ  grain  size  and  age  of  deposition.  The  controlling  factor  in  the  grain 
size  appears  to  be  the  gradation  of  the  fines  (materials  <0.076  mm  particle 
diameter).  More  cementation  is  present  when  the  fines  are  silty  in  nature,  and 
less  cementation  is  present  when  the  sands  are  clean  (i.e.,  very  little  silt). 
This  observation  is  consistent  with  other  areas  and  is  probably  due  to  the 
lower  permeability  of  the  silty  sands  causing  slower  ground  water  movement 
allowing  more  deposition  of  cement  of  the  numerous  grain  contacts. 

The  effect  of  age  of  deposition  is  fairly  obvious--the  younger  the  material, 
the  less  cementation  present.  However,  these  two  phenomena  are  interactive  and 
are  likely  to  be  more  complex  than  the  simplistic  discussion  would  indicate. 

The  mechanics  of  soil  cementation  is  not  yet  clearly  understood,  and  the 
subject  is  still  under  investigation. 

The  seismic  velocity  of  the  granite,  as  determined  from  acceleration  TOAs, 
is  approximately  9300  feet/second.  This  value  is  in  general  agreement  with  the 
visual  classification  of  the  granite  as  highly  jointed  and  fractured.  This 
value  for  seismic  velocity  is  also  consistent  with  the  results  of  the  seismic 
refraction  work  conducted  on  the  HAVE  HOST  Test  Site. 

A  plot  of  the  time  to  peak  velocity  versus  distance  (Figure  3.7)  shows  that 
the  amount  of  cementation  has  little  effect  on  the  propagation  velocity  of  the 
peak  particle  velocity  (~  800  feet/second).  The  18S  and  31.0  feet  18D  loading 
wave  velocities  are  identical;  the  15.0  feet  18D  velocities  are  noticeably 
slower.  This  reinforces  the  observation  from  Figures  3.5  and  3.6  that  the 
maximum  material  response  was  not  a  function  of  the  amount  of  cementation  present 
in  these  alluvial  materials.  Preliminary  grain  size  analysis  (Ref.  5)  indicates 
that  the  materials  grade  finer  with  depth.  The  31.0  feet  (18D)  gauges  were 
placed  in  25  percent  more  fines  that  the  12.0  feet  (18S)  gauges.  At  this  time, 
the  effect  of  the  grain  size  distribution  on  material  response  is  not  evident. 
This  situation  should  clear  up  when  more  laboratory  test  data  become  available. 
Further  analysis  of  the  TOA  and  time  to  peak  velocity  data  indicates  the  ratio 
of  propagation  velocity  to  seismic  velocity,  commonly  assumed  to  be  1/2  as  an 
initial  guess,  does  not  hold  in  this  case.  A  ratio  of  1/3  to  1/4  would  be  much 
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more  accurate.  This  difference  can  noticeably  affect  ground  motion  predictions 
made  for  this  material. 
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SECTION  4 

MATERIAL  MODEL  DEVELOPMENT 

The  primary  purpose  of  the  CIST  18  experiment  was  to  measure  the  dynamic 
response  of  in-situ  materials  for  the  development  of  parameters  in  material 
models  that  could  be  used  in  HAVE  HOST  and  MAP  calculational  efforts.  These 
models  will  serve  as  inputs  for  the  ground  shock  and  structural  dynamics  soil 
island  calculations  that  will  accompany  the  HAVE  HOST  field  test. 

We  cannot  overemphasize  that  material  modeling  is  a  continuing  process,  and 
the  model  presented  in  this  report  is  one  that  is  available  at  this  time  (January 
1977).  This  model  is  by  no  means  intended  to  be  tlie  final  HAVE  HOST  model. 

The  model  parameters  were  developed  using  AFTON  (a  2-D  finite  difference 
code)  and  WONDY  IV  (a  1-D  finite  difference  code).  AFTON  was  used  principally 
because  of  greater  familiarity  with  the  code  and  more  flexibility  of  its  use. 
WONDY  IV  was  used  in  the  parametric  studies  that  were  performed  to  determine 
the  specific  effect  of  individual  parameter  variations. 

The  discussion  of  the  material  model  is  divided  into  the  hydrostatic 
component  and  the  deviatoric  models. 

4.1  HYDROSTATIC  MODEL 

The  hydrostatic  model  contains  both  density  and  energy  components.  First 
the  density  components  are  described  and  then  the  energy-dependent  terms  are 
added. 

The  pressure-density  relation,  Figure  4.1,  is  described  as  the  function 

Ph  -  f  (y.y*)  (1) 

where  p  =  ^~^o,  p  is  the  initial  material  density,  and  y*  is  the  minimum  of 
Po 

Pa  (a  material  parameter)  and  the  maximum  p  experienced  by  the  material.  For  p 
greater  than  pa,  the  hydrostatic  pressure  is  described  by  the  relation 

P^  =  P3+K^(p-P3)  -  (K^-K^)ps[l-Ag]  (2) 

where  A^  =  exp  [-(p-P3)/psl , 
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P3  is  the  hydrostatic  pressure  at  03, 


K 


z 


C 

z 


(1-v,) 


(3) 


and  K^,  C^,  and  are  material  parameters.  In  particular,  is  the 

unloading  Poisson's  ratio.  For  the  condition  p  less  than  the  hydrostatic- 
pressure  relations  are  divided  into  loading  and  unloading  relations. 

The  loading  relations  are  used  when  y  is  greater  than  or  equal  to  y*  and 
are  defined  by  three  relations.  For  y  £  y^,  the  relation  is 


where 


P,,  =  K  y 
1-1  u 


1/  P  l+V 

K  =  C  •  ^0  •  u 

u  T  — 


u 


Pi  =  Pi/K 


(4) 

(5) 


with  P^  and  material  parameters.  In  particular  is  the  seismic  velocity. 
For  yj<y£y2  the  relation  is 


Ph  =  Pi  +  -  Pj)  (6) 

where 

K.  =  C,  •  PQ  •  (7) 

T 

P2  =  Pi  (P2  "  Pi)/\i 

with  Pj,  C|^,,  and  material  parameters.  In  particular,  is  the  loading 
Poisson's  ratio.  For  P2  <  P  £  P,*  the  relation  is 

P^  =  P2  +  -  vj  (8) 

where  ’  f  3:^ 
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with  C^2  ^  material  parameter. 

The  unloading  relation  is  used  when  y  is  less  than  y*.  The  relation  is 
defined  as  a  function  of  the  hydrostatic  pressure,  calculated  by  substituting 
y*  for  y  in  relations  (4),  (6),  or  (8),  and  the  value  of  y.  The  relation  is 


P  =  P  - 

for 

Pp  1  P  <  P 

where 

\ 

For 

M  <  P„ 

where  is  the  local 

lithostatic  load 

(9) 


(10) 


K„  (M-M,) 


2  ■  1 


Two  energy-dependent  expressions  are  used,  depending  on  the  relation  of 
y  to  y^.  For  y  ^  y^,  the  relation  is 


P  =  +  Tpe 


(11) 


B 

where  r  =  A  +  e  +1 

2 

^  % 

P 

n  -p- 
0 

e  is  the  specific  internal  energy  and  A,  B,  and  e^  are  material  parameters. 

c 

For  y  <  y^,  the  relation  is 

P  =  PY'e'  +  Pj^  (12) 

p  =  D^(y  +  1) 
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n 


r 


p/p 


r 


V'  =  A  +  (r-A)(r,^)‘/' 


e'  =  e  -  [1  -  ex|:  ((}))] 

<j)  =  M  (1  -  l/n^)/n^ 

K 


witi)  e,  a  naterial  parameter. 

4.2  DEVIATORIC  STRESS  MODEL 

The  Jeviatoric-stress  relations  are  calculated  using  an  elastic-plastic 
analysis  with  a  nonassociative  flow  rule  for  alluvium  and  an  associative  flow 
rule  for  rock.  The  shear  modulus,  G,  is  calculated  by  the  relation 


where 


6  =  min  (K*A, 


m 


l-2v 

u 

l^^u 


K  is  the  appropriate  K  for  the  loading  condition  and 


•> 

1-2V|_ 

for  p  ^  p 

2 

1+Vj^ 

3 

l-2v 

u 

for  p  <  p 

2 

1+v 

(13) 


The  yield  conditions,  Figure  4.2,  are  functions  of  both  total  hydrostatic 
pressure,  with  compressive-stress  positive,  and  specific  energy.  The  two 
yield  relations  are  expressed  as  follows: 
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[o,C.  (1  -  -^)]  for  P-j.  <  0  (14) 

[C.+S.P^.,  Y,]  forP^>0 

and  i  =  1  is  the  initial  yield-condition,  i  =  2  is  the  residual  yield-condition, 
and  T. ,  C.  S.,  and  Y.  are  material  parameters.  The  nonenergy-dependent  yield 

1  1  j  1  1 

relation  for  any  zone  is  then  calculated  with 

Y'  =  Y'  +  (Y-  -  Y').  (15) 

where  F^  =  min  (0.04  N^,  1) 

and  is  the  number  of  calculational  cycles  of  plastic  flow  for  each  zone. 

Also,  if  zero  is  ever  calculated  from  the  use  of  relations  (14)  and  (15)  or 
if  e  is  ever  greater  than  e^,  then  F^  is  set  to  1.  The  energy  dependence  of  the 
yield  condition  is  then  calculated  with  the  expression 

Y  =  max  (0,  Y'-(l  -  f~)]  (16) 

s 

The  material  parameters  for  the  CIST-18S  and  all  depths,  except  the  69-foot 
depth  in  CIST-18D,  are  given  in  Table  4.1.  The  69-foot  depth  was  not  included 
because  the  material  was  the  same  as  the  62-foot  depth.  Since  the  tests  did  not 
extensively  exercise  material  above  3,000  psi,  any  parameters  that  are  important 
only  above  that  stress  level  are  not  based  on  the  CIST  data.  These  parameters 
(which  include  C  ,  y  ,  K  ,  A,  B,  e  ,  e  ,  Y.,  and  Y^)  were  assigned  values 

based  on  (1)  a  desire  to  be  generic  in  nature,  (2)  previous  experience  that  they 
are  consistent  with  high  pressure  data,  and  (3)  consistency  with  certain  simple 
relations.  Where  those  values  are  important,  extreme  caution  should  be  used, 
and  other  sources  must  be  consulted.  Also,  even  for  the  parameters  within  the 
stress  range  of  interest,  final  agreement  with  CIST  data  (Appendix  B)  has  not 
been  reached.  Therefore,  changes  in  those  parameters  or  material  model  forms 
may  occur. 


max 


mm 


where  P.^.  =  P  + 
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Table  4.1  HAVE  HOST  material  models 


Parameter 

Uni  ts 

18S 

18D 

Depth 

feet 

12 

15 

31 

62  -  ( 

Mat  type 

al  luvi  urn 

al luvi urn 

alluvium 

Rock 

Po 

gm/cc 

1.85 

1.85 

1.85 

2.5 

VL 

0.3 

0.3 

0.3 

0.23 

0.3 

0.3 

0.3 

0.23 

'll 

feet/second 

1500 

1400 

2000 

8500 

C|_2 

feet/second 

1200 

1000 

1400 

8500 

feet/second 

1800 

1800 

2000 

8500 

feet/second 

1800 

1800 

2000 

8500 

psi 

18 

0 

0 

10.00 

psi 

200 

125.00 

125.00 

100.00 

0.17 

0.17 

0.17 

0.04 

0.25 

0.25 

0.25 

0.25 

psi 

E7 

E7 

E7 

E7 

A 

0.5 

0.5 

0.5 

0.5 

B 

1.3 

1.3 

1.3 

1.3 

"o 

cgs 

E9 

E9 

E9 

E9 

cgs 

Ell 

Sll 

Ell 

Ell 

Ti 

psi 

0 

-10.00 

-10.00 

-5.00 

psi 

0 

10.00 

10.00 

5.00 

Si 

0.6 

0.6 

0.6 

1.00 

Yl 

psi 

4000 

4000 

4000 

15000 

^2 

psi 

0 

0 

0 

0 

^2 

psi 

0 

0 

0 

0 

^2 

psi 

0.6 

0.6 

0.6 

1.00 

V2 

psi 

4000 

4000 

4000 

15000 
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Figures  4.3  through  4.6  are  comparisons  of  the  3,  5,  and  8  feet  experiinen- 
tal  data  and  AFTON  calculations  using  the  material  models  given  in  Table  4.1. 

A  similar  comparison  for  12  feet  range  is  shown  in  Figure  4.7.  More  detailed 
comparisons  are  found  in  Appendix  B. 

Comparison  of  the  experimental  and  calculated  data  points  out  several 
deficiencies  in  the  model  currently  used  at  AFWL.  This  model  is  known  as  the 
engineering  model.  The  present  AFWL  computer  codes  (AFTON  and  WONDY  IV)  do  not 
replicate  the  smooth  peaks  in  the  velocity  waveforms  for  dry  materials.  (This 
appears  to  be  a  shortcoming  of  the  engineering  model  rather  than  a  shortcoming 
of  the  codes.)  The  engineering  model  works  with  varying  degrees  of  success  for 
different  materials.  In  dry  sandy  alluvium,  the  deficiencies  are  apparent 
(Figure  4.3). 

In  dry  sandy  alluvium,  the  engineering  model  cannot  completely  match  the 
data  at  the  gauge  ranges  in  the  CIST  (3,  5,  and  8  feet).  The  usual  result 
would  be  to  undershoot  the  data  at  3  feet  and  overshoot  the  data  at  8  feet. 
These  deficiencies  may  be  caused  by  strain  rate  effects  present  in  these 
materials  that  cannot  be  taken  into  account  with  the  engineering  model.  At 
this  point  in  time,  it  is  not  apparent  that  any  other  material  model  can 
adequately  model  this  material  either.  This  point  reemphasized  the  need  to 
develop  material  models  that  can  adequately  quantify  material  response  for  the 
wide  variety  of  materials  encountered  in  MX  geologies. 
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Figure  4.3.  Comparison  of  experimental  data  and  AFTON  calculation 
for  CIST  18S. 
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Figure  4.6.  Comparison  of  experimental  data  and  AFTON  calculation 
for  62*  depth  in  CIST  18D. 
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SECTION  5 
CONCLUSIONS 


Listed  below  is  a  summary  of  the  data  recovery  for  CIST  18. 


Alluvium 
Rock 
TOTAL : 


#  Accel  Msmts  (%  Data  Recovery) 
15  (100) 

12  (100) 

27  (100) 


#  Stress  Msmts  \%  Lata  Recovery) 
8  (100) 

None  recorded 
8  (100) 


It  can  be  seen  that  there  was  100  percent  data  recovery.  Based  on  the  percentatje 
and  quality  of  data  recovered,  CIST  18  was  a  completely  successful  experiment. 
However,  the  data,  wnen  compared  to  the  material  model  calculations,  point  out 
apparent  inadequacies  in  the  engineering  model.  This  problem  is  currently 
being  addressed  at  AFWL/DES. 

Based  on  the  data,  several  conclusions  can  be  drawn  from  the  data  analysis:, 

(1)  Seismic  velocity  is  a  function  of  soil  cementation. 

(2)  The  propagation  velocity  of  the  peak  particle  velocity  is 
independent  of  soil  cementation. 

(3)  The  maximum  alluvial  response  to  CIST  type  loadings  is  independent 
of  soil  cementation. 


This  experiment  has  yielded  a  great  deal  of  information  on  the  material 
properties  of  sandy  alluvium.  However,  much  work  remains  to  be  accomplished, 
particularly  in  the  development  of  a  material  model  to  adequately  characterize 
the  response  of  these  materials. 
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APPENDIX  A 

ACCELERATION,  VELOCITY  AND  DISPLACEMENT  TIME  HISTORIES 

Plots  are  arranged  in  order  of  increasing  range  (3,  5,  8,  and  12  feet, 
respectively).  Horizontal  data  for  a  given  location  are  followed  by  vertical 
data  if  available.  Whenever  available,  redundant  recordings  are  given  in 
lieu  of  primary  recordings.  The  redundant  recordings  were  calibrated  at 
80  percent  band  edge,  and  the  primary  recordings  were  calibrated  at 
40  percent  band  edge;  hence,  tne  redundant  recordings  signal  to  noise  ratio 
is  much  improved. 

A.l  MEASUREMENT  IDENTIFICATION 

Each  data  trace  is  identified  at  its  top  center  by  a  measurement  designation 
number.  The  measurement  designation  number  consists  of  eight  alpna-numeric 
designators  in  the  following  form: 

X  -  X  -  XXX  -  XX. X  -  XXX  -  X.X  -  XX  -  X 

1  2  3  4  5  6  7  8 

a.  The  first  character  indicates  the  organization  that  established  the 
measurement  required: 

F  -  AFWL  (Free-Field)  (DES-G) 

b.  The  second  character  denotes  the  method  of  data  acquisition: 

E  -  Electronic 

c.  The  third  set  of  characters  indicates  the  plan  location  of  the 
free-field  measurement.  LOl  refers  to  hole  number  one  (hole  layouts  are 
given  in  Figures  2.6  and  2.7). 

d.  The  fourth  set  of  characters  indicates  the  depth  (in  feet)  of  the 
transducer  below  the  surface. 

e.  The  fifth  designator  indicates  the  azimuth,  in  degrees  from  North 
(0°),  of  the  radial  on  which  the  measurement  is  made. 

f.  The  sixth  set  of  characters  indicates  the  radial  distance  in  feet 
from  GZ  to  the  center  of  the  transducer. 

g.  The  seventh  set  of  characters  specifies  the  type  of  measurement 
being  made  (uS  shown  below). 

A  -  Acceleration 
CP  -  Cavity  Pressure 
SE  -  Soil  Stress 
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h.  The  last  set  of  characters  indicates  the  orientation  of  the  sensing  axis 
of  the  transducer: 

V  -  Vertical 
H  -  Horizontal  Radial 
T  -  Horizontal  Transverse 

Following  the  measurement  designation  number  are  the  tape  number,  track 
number,  and  measurement  number.  These  three  numbers  appear  to  the  right  of  the 
measurement  designation  number  and  are  separated  by  a  slash  (/)  and  an  asterick 
(*).  The  first  is  the  tape  machine  number;  the  second  is  the  track  number,  and 
the  third  is  the  measurement  number. 

Figure  A.  1  illustrates  the  azimuti,  range,  and  sensing  axis  of  the  free- 
field  transducers.  A  sample  of  the  labeling  system  is  given  in  Figure  A. 2. 

Table  2.1  presents  the  CIST  measurement  list  and  calibration  values. 

A. 2  DATA  PLOTS 

The  time  histories  presented  in  the  appendix  are  edited  data.  The  only 
posttest  data  correction  that  has  been  made  is  that  cable  breaks  were  truncated. 
This  was  done  so  the  first  and  secord  integral  time  histories  scaling  would  not 
be  driven  by  the  cable  breaks.  Data  correction  is  planned  in  the  near  future. 

It  appears  at  this  time  that  there  are  no  drastic  early  time  (first  25  msc) 
corrections  to  be  made. 


4  6 


AFWL-TR-77-81 


/ 


G1 


L02 


GAGE  ORIENTATION 


AFWL-TR-77-81 


48 


Figure  A. 2.  Sample  of  record  labeling  system. 
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APPENDIX  B 

MATERIAL  MODEL  COMPARISONS 


C  -  Calculated  Horizontal  Velocity 
V  -  Calculated  Vertical  Velocity 
E  -  Experimental  Horizontal  Velocity  (Uncorrected) 
0  -  Experimental  Vertical  Velocity  (Uncorrected) 
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Figure  B.l.  Comparison  of  experimental  data  and  AFTON  calculation 
for  3'  range  and  12'  depth  in  CIST  18S. 
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Figure  B.2.  Comparison  of  experimental  data  and  AFTON  calculation 
for  5'  range  and  12'  depth  for  CIST  18S. 
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Figure  B.3.  Comparison  of  experimental  data  and  AFTON  calculation 
for  8'  range  and  12'  depth  for  CIST  18S. 
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Figure  B.5.  Comparison  of  experimental  data  and  AFTON  calculation 
for  5'  range  and  15'  depth  for  CIST  18D. 
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Figure  B.6.  Comparison  of  experimental  data  and  AFTON  calculation 
for  8'  range  and  15'  depth  for  CIST  18D. 
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Figure  B.7.  Comparison  of  experimental  data  and  AFTON  calculation 
for  3'  range  and  31'  depth  for  CIST  18D. 
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Figure  B.8.  Comparison  of  experimental  data  and  AFTON  calculation 
for  5'  range  and  31'  depth  for  CIST  18D. 
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Figure  B.9.  Comparison  of  experimental  data  and  AFTON  calculation 
for  8'  range  and  31'  depth  in  CIST  18D. 
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Figure  B.IO.  Comparison  of  ijxperimental  data  and  AFTON  calculation 
for  3'  range  and  62'  depth  for  CIST  18D. 
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Figure  B. 12. 


Comparison  of  experimental  data  and  AFTON  calculation 
for  8'  range  and  62'  depth  for  CIST  18D. 


97 


AFWL-TR-77-81 


Figure  B.13.  Comparison  of  experimental  data  and  AFTON  calculation 
for  12'  range  and  62'  depth. 
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Figure  B.16,  Comparison  of  experimental  data  and  AFTON  calculation 
for  5'  range  and  69'  depth  in  CIS!  18D. 
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Figure  B.17.  Comparison  of  experimental  data  and  AFTON  calculation 
for  8'  range  and  69'  depth  in  CIST  18D. 
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Figure  B.18.  Comparison  of  experimental  data  and  AFTON  calculation 
for  12'  range  and  69'  depth  in  CIST  18D. 
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